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1. Introduction . In our first paper on this subject (Fan and Wu, 1981), we derived 
an expression for the rate of energy generation by a pulsar and estimated the 
contribution from all the pulsars in our galaxy to the observed cosmic ray 
intensity. The theory was then developed to an expanded version, and 
observational facts supporting the theory were cited (Fan et al., 1982; Fan and 
Hang, 1983). In this paper we supplement additional evidences. 


2. Brief Review . The theory is based on two propositions: 

(1) The rate of the energy generation by a spinning pulsar is given by 


— = K M2(1)2 


dt 


Rc 


( 1 ) 


where M, R, and to are respectively the magnetic dipole moment, the radius, and 
the angular velocity of the pulsar; c is the velocity of light, and K is a 
characteristic constant of the surrounding medium. This expression should be 
viewed as Ohm's Law in cosmological scale in which the EMF is generated by the 
spinning magnetized celestial body and K is a measurement of conductivity of the 
medium. 


(2) The magnetic dipole moment of the pulsar, M, is related to its angular 
momentum L by the following expression: 


M = otL , (2) 

where a is a "form factor," a characteristic constant of the interior of the pulsar. 
It is an extrapolation of so-called "Bode's Law" of planetary magnetism. The 
relationship seems to hold approximately true for all celestial bodies whose 
magnetic diple moments and angular momenta have been estimated or measured. 
They are over 15 orders of magnitude with Mercury, Venus, and Mars at the lower 
end and 24 Babcock magnetic stars at the upper end of the ladder, including a 
neutron star in the middle. The physical reason behind this fundamental relation 
is not known. 

By combining (1) and (2), we have for the rate of energy generation as 

dT _ Kg 2 L 2 to 2 _ Ka 2 1 2 to 4 
dt = c 2 R ~ c 2 R ’ 

where I is the moment of inertia of the pulsar. Since T = 1/2 Ico 2 , we have finally 

_ _2_ do) 2Ka 2 I 
u) 3 dt 


c 2 R 


(4) 
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The values of K and a were determined from the observation on Jupiter; they are 
2xl0“ 2 /sec and 2x10" 16 gauss cm.s.g' 1 , respectively. The scaling constant kot 2 is 
then 8x10“ 3< * cm.s" 1 .g* 1 . Equation (4) allowed us to use measured values of w and 
dw/dt of a pulsar to calculate its I/R, and then, by invoking neutron star models, 
we determined the mass m, the moment of inertia I, and the radius R of the 
pulsar and calculated the value of dT/dt. 


3. Additional Evidences . An integration of Eq. (4) yields 

1 1 t 2Kot 2 I 

w 2 - w 2 = Yt , Y - c 2 R , 


(5) 


where a> 0 is the initial angular velocity of the pulsar. Assume that to 0 >>03, then 

\ = yt . ( 6 ) 

or 

Figure 1 is a plot of log(2/w 3 ) (-dw/dt) against logo) for 202 pulsars (Fan et al., 
1982). The pulsars at the left lower corner of the figure appear to be on a 
straight line given by Eq. (6) with t as a parameter. We interpret the line as the 
age limit of all pulsars. Using 8x10“ 34 cm.s^g" 1 for the value of Kot , we find t ~ 
6x10 15 sec (0.5 billion years). Consider the fact that the scaling constant could be 
off by a factor of 10, the time limit seems to be reasonable. 

The low-energy charged particle (LECP) on Voyager 1 and 2 also detected a 
flux of hot plasma escaping from the Saturnian magnetosphere, just like in the 
case of Jupiter (Krimigis et al., 1981). Although the dynamics of the Saturnian 
magnetoshpere may be complicated by the existence of the rings, the energetic 
particles escaping from the system can be used to make an order-of-magnitude 
check of the prediction of Eq. (1). Estimating from their published figures, the 
escaping flux of 53-85 keV ions from the Saturn appears to be a factor of 100 
lower than that from the Jupiter. The prediction of Eq. (1) is 1000. 
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